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We investigate the translocation of a single stranded DNA through a pore which fluctuates between 
two conformations, using coupled master equations. The probability density function of the first 
passage times (FPT) of the translocation process is calculated, displaying a triple, double or mono 
peaked behavior, depending on the interconversion rates between the conformations, the applied 
electric field, and the initial conditions. The cumulative probability function of the FPT, in a 
field-free environment, is shown to have two regimes, characterized by fast and slow timescales. 
An analytical expression for the mean first passage time of the translocation process is derived, and 
provides, in addition to the interconversion rates, an extensive characterization of the translocation 
process. Relationships to experimental observations are discussed. 
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INTRODUCTION 

Translocation of biopolymers through a membrane 
pore occurs in a variety of biologic al processes, such 
as ge ne expression in eucaryotic cells l)B. Alberts et al. L 
Il994j) . co njugation betwee n proc aryotic cells, and virus 
infection l|Madigan et al. L Il997j) . The importance of 
translocation in biological systems and its applications, 
have been the motivation for recent theoretical and ex- 
perimental work on this topic. In experiments one usu- 
ally measures the time it takes a single voltage-driven 
single stranded DNA (ssDNA) to transloca te through a- 
hemo lysin channel of a known structure l|Song et al. L 
1996). Since ssDNA is negatively charged (each monomer 
of length b has an effective charge of zq, where q is the 
electron charge, and z (0 < z < 1 ) is controlled by 
the solution pH and strength), when applying a voltage 
the polymer is subject to a driving force while passing 
through the transmembrane pore part (TPP) from the 
negative (cis) side to the positive (trans) side (for an 
illustration of the process see Meller 2003, figure 4). Be- 
cause the presence of the ssDNA in the TPP blocks the 
cross- TPP current, one can deduce the FPT probability 
density functio n (pdf), F(i), from the current blockade 
durat ion times ijKasianowicz "eTlLTIilMfMeller et al. L 
1200 ll) . 

Experiments by Kasianowicz et al. (1996), show F(t) 
with three peaks. It was suggested that the short- 
time peak represents the non-translocated events, while 
the other two, longer-time peaks, represent transloca- 
tion events of different ssDNA orientations. In addition, 
the times that maximize the translocation peaks were 
shown to be proportional to the polymer length and in- 
versely proportional to the applied field. In experiments 
by Meller et al. (2001), F(i) was shown to be mono- 
peaked, with a corresponding maximizing time that has 
an inverse quadratic field dependence. More recently, 



Bates et al. (2003) measured the FPT cumulative prob- 
ability function (cpf), which is the probability to exit the 
channel until time t, G(t) = J Q F(s)ds, in a field-free en- 
vironment. G(t) was approximated by two well separated 
timescales with the ratio of 1/20. 

In previous theoretical works, the translocation of a 
ssDNA through a nanopore was described by statisti- 
cal models that focused on calculating the free energy of 
the process as a function of the translocation state. The 
free energy contained terms that represent the entropy 
and the chemical potential of the polymer parts on both 
sides of a zero thickness membrane (Sung and park, 1996; 
Muthukumar, 1999). The role of the membrane thickness 
was studied by Ambjornsson et al. (2002), Solonika and 
Kolomeisky (2003), and Flomenbom and Klafter (2003a). 
The obtained free energy was used to calculate the mean 
first passage time (MFPT), which asymptotically was 
found to scale linearly with the polymer length for a field- 
biased process. This is the expected MFPT dependence 
of a Markovian biased random walk in a finite interval 
(iRednerLElol. 

A different approach was suggested by Lubensky and 
Nelson (1999) and was further developed by Brezhkovskii 
and Gopich (2003), where a diffusion-convection equa- 
tion was used to describe the translocation process, un- 
der the assumption that the polymer parts outside the 
membrane hardly affect the translocation. Brezhkovskii 
and Gopich (2003) showed that by changing the cis ab- 
sorbing end to be partially absorbing, the mono-peaked 
F(t) obtained by Lubensky and Nelson (1999) changes to 
a superposition of a decaying non-translocation pdf and 
a peaked translocation pdf. Using the fractional Fokker- 
Planck equation, Metzler and Klafter (2003) suggested 
an explanation for the slow relaxation time of the exper- 
imentally obtained G(t). We have shown by using the 
master equation (ME) that F(t) can be double or mono 
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peaked, de pending; on the applied field and on the initial 
condition ijFlomenbom and Klafterl l2003a|) . 

In the approaches summarized above the structure of 
the pore is taken to be rigid, namely, governed by a 
single conformation. Although it is known that the a- 
hemolysin cha nnel has a rigid st ructure that allows its 
crystallization ijSong et al. Ul996|) . during the transloca- 
tion of a long polymer (larger than the pore length) with 
the same width that of the channel at some cross sec- 
tions along the channel, small fluctuation in the channel 
structure may occur, which give rise to a more complex 
process than what was assumed so far. In this work we 
relax the assumption of a single pore conformation and 
introduce a second conformation coupled to the first one. 
In a continuum formalism, the process takes place in 
an effective two dimensional system, where one dimen- 
sion represents the translocation itself, and the second 
dimension represents the structural fluctuations. This 
picture is richer and is more realistic, since small struc- 
tural changes in physiological conditions are known to 
occur in large biomolecules, certainly during interaction 
with other large biomolecules. 

The function that represents best the translocation 
process is F(t) (or its integral G(t)). Through the depen- 
dence of F(t) on the system parameters, we learn about 
the important degrees of freedom which participates in 
the translocation process. The characteristics of F(t) are 
the dependence of its shape, moments, and times that 
maximize its peaks on the system parameters. Using the 
generalized model that takes into account fluctuations in 
the pore structure, we calculate F(t) and show that it can 
display one two or three peaks, depending on the applied 
voltage, the temperature and the interconversion rates 
between the two conformations. Analytical expressions 
for the MFPT are derived and related to the experimental 
findings. In addition, we calculate the cumulative prob- 
ability G(t) in the field-free limit, and show that it also 
provides valuable information about the system param- 
eters. Thus, these tools help in gaining insight into the 
translocation of a polymer through a narrow pore, and in 
explaining the diversity of the experimental obs ervations 
l|Kasianowicz et al. I Il99l iMeller et al. Ll200l|) . 

THEORETICAL MODELLING 

Basic Model 

The basic model we use to describe the translocation 
relies on a one-dimensional process. To use this sim- 
plification, we map the three-dimensional translocation 
process onto a discrete one-dimensional space containing 
n(= N + d— 1) states separated from each other by a unit 
length b. The translocation takes place within a TPP of 
a length that corresponds to d monomers. An n-state 
ME is introduced to describe the translocation of an N- 
monomer long ssDNA subject to an external voltage V 
and temperature T. The occupation pdf of the j state 
is [P (t))j=Pj(t), where the state index j determines the 
number of monomers on each side of the membrane and 



within the TPP (rrij ) . Pj (t) satisfies the equation of mo- 
tion: 

dPj{t)/dt = a j+1>j P j+1 (t) + Oj-xjPj-xit) - 

under absorbing boundary conditions on both sides of 
the membrane (the polymer can exit the TPP on both 
sides). Eq. 1 can be written in a matrix representation: 

dP(t)/dt = A~P(t), (2) 

where the propagation matrix A is a tridiagonal matrix 
that contains information about the transitions between 
states in terms of rate constants, a 3 ; j±i, which are given 
by: 

: I = / '../';„ : I- ( 3 ) 

Here kj is the rate to perform a step, Pjj-i (pjj+i) is the 
probability to move one state from state j to the trans 
(cis) side, and Pjj+i +Pjj-\ = 1. 

kj is taken to be similar to the longest bul k relaxation 
time of a polymer ijDoi and Edwards! 1 1986|) : 

k 3 ; = l/(/3^6 2 mp = R/m^ = k B T, (4) 

with two exceptions: the parameter £ p represents the 
ssDNA-TPP interaction and cannot be calculated from 
the Stokes relation, and \x serves as a measure of 
the polymer stiffness inside the confined volume of 
the TPP, and is bounde d by the conventional values 
i|Doi and Edwards I. Il986|) : < /x < 1.5. 

Assuming a quasi-equilibrium process, which enables 
using the detailed balance condition, and using the ap- 
proximation cij.j-i/ 'aj_ij ~ pjj—i), the prob- 
ability Pj.j-i is found to be: 

p j j-i = 0- + e^)- 1 . (5) 

The free energy difference between states, AEj = Ej_i — 
Ej, is computed considering three contributions: electro- 
static, entropic, and an average attractive interaction en- 
ergy between the ssDNA and the pore. More explicitly, 
f3AEj is given by 0AEj = 0AE? + Sj, where /3A£f < 
represents the effect of the field which directs towards the 
trans-side and 6j > (for j > d) represents an effective 
directionality to the cis-side, which originates from the 
entropic factors and the average attractive interaction 
energy between the ssDNA and the pore. For a more 
detailed discussion see Flomenbom and Klaftcr (2003a). 

Several features emerge from the simple one- 
dimensional model. For homopolymers, poly-dnu, where 
nu stands for the nucleotide type, we estimate £ p (A nu ) ~ 
10~ 4 77ieUs/?im 2 , £,p{C nu )=£ p {T nu )=£, p {A nu ) /3 and 
H(C nu )=l, n(A nu )=l.U, fi(T nu )=1.28. Here A nu , 
C nu and T nu stand for adenine, cytosine and thymine 
nucleotides, respectively. Interestingly, £ p is three order 
of magnitude larger than the bulk friction constant, 
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which is consistent with the role assign to £ p to represent 
the interaction between the polymer and the channel. 

In addition, from the expressions for pAEj and pjj^i, 
the important parameter V/Vc = j3z\q\V(l + l/d) comes 
out naturally. This ratio determines the directionality of 
the translocation, and, in particular, for V/Vc > 1 there 
is a bias towards the trans-side of the membrane. 

Translocation Through a Conformationally 
Changing Pore 

A more realistic description of the translocation can 
be obtained by taking into consideration fluctuations in 
the TPP, either spontaneous or interaction induced. Ac- 
cordingly, we introduce an additional pore conformation 
which is represented by the propagation matrix B. The 
changes in the pore conformation between A and B are 
controlled by the interconversion rates, uia and ujb- 
(lob) is the rate of the change from the A (B) to the B 
(A) pore conformation. 

The physical picture of the process is that when the 
pore conformation changes, a different environment is 
created for the ssDNA occupying the TPP. This implies a 
change in £ p and /i. For a large polymer, N > d, we take 
B w AA, where A is a (dimensionless) parameter that 
represents the effect of the conformational change on £ p 
and /j,. The parameter A may be interpreted as a measure 
of an effective available volume in the TPP, when the 
amino acids residues protruding the TPP change their 
positions. 

The equations of motion of the ssDNA translocation 
through the fluctuating pore, written in matrix represen- 
tation, are: 

d_ (p(t;A)\ = (A-uiA c B \(p(t;A)\ 

dtyp(f,B)j v Wa b-^b J y ~p(t-B) ) ' 

(6) 

where P(t;i),i = A, B is the occupation pdf vector of 
configuration i, oj\ = u>jl, and I is the unit matrix of n 
dimensions. 

As a general note we refer to the form of Eq. 6, 
which was used to study the resonant a ctivation phe- 
nomenon (B ar-Haim an d Klaftcr 1 119991) . This phe- 
nomenon, which was first reported by Doering and 
Gadoua (1992), is the occurrence of a global minimum 
in the MFPT as a function of the interconversion rate 
for a system in which lua = u)b- Because of the assump- 
tion B = AA, the s ystem investigated here cannot e xhibit 
this phenomenon l|Flomenbom and Klafterl |2004|) . 

DENSITY OF TRANSLOCATION TIMES 

Parameters Tuning 

To study the translocation of ssDNA through a fluctuat- 
ing pore, we start by computing F(t). Formally, F(t) is 
defined by 

F(t) = 9(1 S(t))/dt. (7) 



Here, S(t) is the survival probability, namely, the proba- 
bility to still have at least one monomer in the pore, and 
is given by summing the elements of the vector that solves 
Eq. 6; see Appendix A for details. Using the known val- 
ues of £ p and [i from the single conformation model, we 
examine in this subsection the effect of the parameters 
A, uja and lub on F(t). 

First, we check the effect of A on F(t) for several lim- 
iting cases. For A=0 movement in any direction occurs 
only under the A conformation environment. The B con- 
formation traps the polymer for a period of time governed 
by the interconversion rates. For A=l, namely, B=A, the 
environmental changes do not affect the translocation, 
and the process reduces to a translocation through a sin- 
gle conformation. For A > 1 the environmental changes 
enhance the process. In this paper we restrict ourselves 
to the range < A < 1. 

The picture is less intuitive for intermediate values of 
A. Fig. 1 shows that by choosing A properly, three peaks 
in F(t) can be obtained. In particular, as shown in the 
inset of Fig. 1, the range of A values for which F(t) 
exhibits three distinct peaks is 0.1 < A < 0.30. For 
the single conformation case we found that F(t) can be 
either mono or double peaked depending on V/Vc, an d 
on the initial state of the translocation x. The short 
time peak represents the non translocated events, while 
the long time peak represents the translocation events. 
The generalization for two pore conformations may yield 
two translocation peaks in addition to a short time non 
translocation peak. Indeed, Fig. 1 supports the expected 
behavior for the limiting A values, and shows that as 
A — ► 1, F(t) possess only one translocation peak, as well 
as for A — > 0, where the B conformation peak spreads out 
towards larger times, which results in its disappearance. 

Although Fig. 1 is obtained for a given value of the in- 
terconversion rates, our explanations regarding the F(t) 



,1 10 

Si — i — 




0.3 1 ti 2 2.5 3 

t(ms) 



Figure 1 Poly-dT ntl F(t), for several values of A, with: JV=30, 
d=12, x=N + d/2, T=2°C, V/V c =2, w s =10 2 Hz, w=l, and 
z ~ 1/2. The left peak represents the non-translocated 
events, whereas the other two peaks represent translocation. 
Inset: The range for which A yields three peaked Fit) is shown 
to be 0.1 < A < 0.3, when given the above parameters. 
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behavior for the limiting cases A=1,0, are valid for any 
system conditions. This is demonstrated by calculating 
the MFPT (Appendices B and D). In Appendix B we 
show that when A=l, the MFPT of the two configura- 
tions model reduces to that of the single conformation 
model. In Appendix D we show that for A=0, the B con- 
formation contribution for the MFPT is a term which is 
inversely proportional to the interconversion rate, lu^ 1 . 

Therefore, A serves as a tuning parameter that leads 
to either one or two actual translocation peaks in F(t). 
The question of interest is how A depends on the system 
parameters. We assume a small field perturbation in the 
regime of biological interest [0 < V/Vg < 3, using Vg re 
50mV l|Flomenbom and Klafterl . l2003atv j . so that X(V) 
follows A re Ao + V/V\, and keeping X(V) < 1. Here Ao 
and V\ might be expansion coeffi cients, where An «C 1 is 
implied from recent experiments ((Bates et al. 1 12003^ . as 
we discuss later. The process can be viewed such that 
as the voltage increases, amino acids residues protruding 
the TPP that constitute obstacles for the translocating 
ssDNA clear the way. While the A dependence on the 
voltage is assumed here, its dependence on other system 
parameters (e. g., temperature and pH) is unknown and 
is folded into V\. 

To check how interconversion rates affect F(t), it is 
convenient to define two dimensionless parameters, lu = 
to a/ to b and tus/k (or LUA/k), where k is the dominate 
rate of the A conformation for a sufficiently large N, 
k = R/d* 1 . The first ratio set the dominance of a given 
conformation over its counterpart; namely, for lu <C 1 
most of the translocation events take place in the A con- 
formation. The second ratio gives an estimate of the 
number of moves done in a given conformation before a 
change in the pore structure occurs, and thus relates the 
ssDNA dynamics to the structural changes dynamics. 

As shown in Fig. 2 and in the inset of Fig. 2, F(t) 
exhibits two peaks corresponding to actual translocation 
only when u> re 1. For lu <C 1 and lu 3> 1 only one peak 
corresponding to an actual translocation survives. For 
all cases there is a peak representing non-translocation 
events. In addition, we find that for two translocation 
peaks to be obtained, the ratio lob / k (or uiA/k due to lu re 
1) must fulfil LUs/k < 1CP 3 (data not shown). The lower 
limit of the interconversion rates is inversely proportional 
to the order of the measurement time otherwise only one 
conformation will be detected. 

Finally, we assume that the rate of the conformational 
changes is controlled mainly by temperature; namely, we 
take lua and wb as voltage independent in the regime of 
biological interest: < V/Vg < 3. 

Translocation Velocity 

To study further the translocation process, we check 
the voltage dependence of the times that maximize the 
peaks of F(t), denoted as t m ^ where i=l, 2, 3 (e.g. t m A 
characterized the short time pea k ). In previou s work s 
ijFlomenbom and Klafter 1 l2003at iMeller et al. j I2001D . 
f" 1 for one translocation peak was regarded as the most 
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Figure 2 Po\y-dT nu F(t), for several values of lua and fixed lub 
(o;_b = 10 2 Hz), with A = 1/4, and the other parameters as in 
Fig. 1. Inset: For small values of lu, lu < 1CT 2 , F(t) displays 
one translocation peak that corresponds to conformation A, 
whereas for large values of lu, lu > 10 2 , F(t) displays one 
translocation peak that corresponds to conformation B. For 
lu ~ 1 two translocation peaks are obtained. 

probable average velocity of the translocation (up to a 
multiplicative constant). We show below that our as- 
sumptions regarding the voltage dependence of the sys- 
tem parameters, yield either linear or quadratic depen- 
dence of the translocation velocity on the voltage, and 
can be used to explain the different experimental obser- 
vations. 

Fig. 3 shows t m 4Vg/V) and t^V/Vc), for 
V\=350mV, in a voltage window that leads to 0.215 < 
A(V) < 0.30, and accordingly to a triple-peaked F(t). 
t m A is almost independent of Vg/V (see Fig. 3a). Al- 
though the non translocation peak amplitude decreases 
upon increasing V/Vg, the location of its maximum 
hardly changes. This happens since exiting against 
the field must happen within a short time window at 
the beginning of the process, otherwise the polymer 
is more likely to cross the membrane due to the elec- 
tric bias. Similar behavior was observed experimentally 
ijKasianowicz et al~T 11996^1 . 

For the single conformation case, we showed that 
t ~ 2 scales linearly with V/Vg when the initial state 
of the translocation is near the cis-si de of the mem- 
brane ijFlomenbom and Klafterl l2003a|) . Fig. 3b shows 
that the linear scaling of t^ 2 (V/Vg) persists. How- 
ever, t^ 3 {V/Vg) (Fig. 3c) displays a quadratic behavior, 
which is a consequence of the form of \{V), as discuss 
in the next section when computing the MFPT. On the 
other hand, setting V\—120mV leads to one transloca- 
tion peak, and to small deviations from linearity towards 
a weak quadratic behavior of ^^(W^c); see Fig. 4. 

The model of two conformations not only yields one or 
two actual translocation peaks as a function of V\, but 
can account for either a linear or quadratic dependence 
of the translocation velocity with the voltage, again, as 
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Figure 3 a: t m> i for poly-dT nu , as a function of V/Vc for 
the same parameters as in Fig. 1 and V\ = 350mV. t m ,i is 
almost independent of Vc/V in contrast to the pronounced 
dependence of t m ,2 and t m ,3- b-c: t m * 2 and t~* 3 depend lin- 
early and quadratically on V/Vc , respectively. The solid lines 
through the circles are polynomial fits. 
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Figure 4 t m 2 for poly-dT mu for the same parameters as in 
Fig. 3 except for V\ = 120mV '. This value for V\ leads to 
0.625 < \(V) < 0.875 and accordingly for one translocation 
peak. The solid line is a polynomial fit. Inset: t m ,i and t m ,2 
behave qualitatively the same as for the case V\ = 350ml 7 . 



a function of V\. Thus, varying V\ we obtain different 
behaviors of the translocation, which can be related to 
the different experimental observations. 

THE MFPT 

Small Field Biased Translocation 

We now turn to calculate the MFPT, which allows for 
an analytical estimation of the characteristic times of the 



FPT pdf and cdf. In general, the m moment of F(t) is 
calculated by raising to the m power the inverse of the 
propagation matrix. For the two conformations translo- 
cation this matrix is given on the right hand side of Eq. 
6. 

After somewhat lengthly calculations, which are given 
is Appendices B and C, the expression for the MFPT, 
< t >, reads: 



< T >r 



A 



[{XPa,o+Pb,o)+t(u a +ub)/2], (8) 



where r is the MFPT for the single configuration model, 
and is given by Eq. C5, and Pa,o(Pb,o) 1S the probability 
that the process starts in conformation A (B). For Pa,o 
and Pb,o the equilibrium condition is assumed, Pa ft = 
lu b / {uJa + w B ) and P Bfi = 1 - P A ,o- 

Eq. 8 is valid for not too high fields, V/Vc ^ 1, and 
the relations between the interconversion with k found 
in the previous section, u>A/k,uJb/k <C 1. The first term 
in the brackets of Eq. 8, XPa,o + Pb,o, represents the 
translocation peaks and can be compared with t m p. and 
tm.3- The second term in the brackets, ~t{loa + Wb)/2, 
represents the coupling time cost, and is of the order of 
o(10~ 2 ) for voltages that obey V/Vc > 1-5. Keeping the 
first term in Eq. 8, we have 



< T >r 



2x£, p b 2 d» 



1 



z\q\(l + l/d)V- 



y-(PA,o + P B ,oy), (9) 



where x « N for a translocation process that starts near 
the cis-side of the membrane. 

Eq. 9 provides a solid basis for the numerically ob- 
tained dependence of the translocation velocity on the 
voltage. < t > consists of two terms that can be at- 
tributed to the A (first term in the brackets) and B 
(second term in the brackets) pore conformations. For 
V\ w 120m V we can replace the expression in brackets 
by one in the relevant voltages window. Thus, we find 
that < t >oc (V— Vc) -1 , which implies that F(t) has one 
translocation peak for this choice of V\. For higher values 
of V\ and voltages of biological interest, the two terms in 
the brackets are separated. This leads to a term that rep- 
resent the A conformation and scales as (V — Vc)~ x , and 
a term that represents the B conformation that scales as 

mv-vc)}- 1 . 

Accordingly, Eq. 9 captures the physical essence of 
the translocation of the ssDNA through the conforma- 
tionally changing pore, under a relatively small field. 



Field-Free Translocation 

In recent field- free experiments by Bates et al. (2003), 
the cdf G{t)=$l F(s)ds was shown to have two regimes 
that were approximated by a fast and a slow timescales, 
n and r 2 , with the ratio ti/t 2 w 1/20. Motivated by 
these experimental result, which implies , within our ap- 
proach, that Aq fulfils Aq <C 1, we study in this subsection 
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the zero field translocation, V — > 0. We start by com- 
puting G(t) for a translocation process that starts at the 
middle state, x = n/2. This is the the sam e initial condi- 
tion t hat was imposed in the experiments ijBates et al. I 
2003). As shown in Fig. 5 (full curve), G(t) displays two 
regimes, a fast increase at short times and a slow increase 
from intermediate to large times. Accordingly, we try the 
approximation 

G ap (t) « 1 - (PA,oe- t/Tl + P Bfi e- t/T2 ). (10) 

Matching the first and the second moments obtained 
from F(t) and from the approximated F(t) we find that 
the characteristics timescales of G ap (t) are (see Appendix 
D): 

ri =t(1 + 3cj/2); t 2 = t{1/2 + lo) + 1/uj b , (11) 

which, when used in G ap (i), lead to the dashed curve 
plotted in Fig. 5. Also shown, dotted curve, is a modified 
version of G ap (t), where t\ — » t m is used in Eq. 10. Note 
that for the short times, t < ti, the later approximation 
fits better G(t), but from intermediate times, t > 3ti, 
G ap (t) and G(t) coincide. 

The two conformations model produces a temporal be- 
havior that agrees with experimental observation, and 
provides a good explanation for it. In the limit, V — > 0, 
the B conformation acts as a trapping conformation; 
namely, the polymer is stuck in its position when subject 
to the environment due to the B conformation. Move- 
ment occurs only through the A conformation. As a re- 
sult two regimes are obtained for G(t). The fast increase 
in G(t) at short times, is a consequence of exiting due 
to the A conformation (at either side of the membrane) , 
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Figure 5 G(t) for poly-cL4 ntl , full curve, for V=0, and the 
initial state x—n/2, ua=10~ 2 'Hz, w=1/2 and the other pa- 
rameters as in Fig. 1. Also shown is the approximate cdf 
G a p(t), dashed curve, and its modified version, dotted curve. 
Inset: F(i) for po\y-dA nu for the corresponding cdf shown in 
the main Fig. 



while the slow saturation at longer times is a result of 
the release from the trapping in the B conformation. 

In the inset of Fig. 5 we show both F(t), the ap- 
proximate F(t), and the modified version of the approx- 
imation, which is obtained when using n — + t m . Be- 
cause the process starts in the middle state, x=n/2, F(t) 
has only one peak, which coincid es with previous results 
IjFlomenbom and Klafter ll2003a|) . Although the approx- 
imate F(t), or any other approximation of two exponen- 
tials with positive coefficients, does not exhibit a peaked 
shape, information about the maximal peak value of F(t) 
and the interconversion rates can still be extracted from 
G ap (t) timescales by using Eq. 11. For example, the 
timescales suggested by Bates et al. (2003), imply that 
t m w 165/iS, and ojb ~ 300Hz. 

CONCLUSIONS 

The model introduced here describes the translocation 
of ssDNA through a fluctuating pore structure. As a 
consequence the ssDNA within the transmembrane pore 
part is exposed to a changing environment, which can be 
reflected in the first passage times pdf, F(t). By comput- 
ing F(t), comparing our results to experimental results, 
and using physical arguments, we obtained theoretically 
a behavior which was previously observed experimentally 
- F(t) having three peaks. This behavior is obtained by 
tuning the dimensionless parameter A, which controls the 
effect of the change in the pore structure on the translo- 
cating ssDNA, and the interconversion rates between the 
pore conformations, to a and ojb- In particular, A has 
to fulfill 0.1 < A < 0.30, and the interconversions rates 
have to be of the same order of magnitude, and much 
smaller than the typical rate of the A pore conforma- 
tion, k, oJB/k < 10~ 3 . This implies that the relaxation 
timescale of the ssDNA in the pore is much shorter than 
of the pore conformational changes timescale. From these 
relations the maximal values of the interconversion rates 
can be deduced from the value of k, given by Eq. 3, to 
be: oja ~ ujb—10 2 Hz. 

We have been able to show, both numerically and an- 
alytically, that the times that maximize the F(t) actual 
translocation peaks, t m ,i,i = 2,3, and the MFPT, are 
inversely proportional to the first or the second power of 
the field, depending on V\. This emphasizes the crucial 
role of V\ on the translocation extracted functions, and 
may explain the different experimental results for F(t) 
discussed in the introduction, meaning that V\ is sensi- 
tive for the specific experimental set up, and biological 
conditions. 

The probability to exit the channel until time t, G(t), 
in a field-free environment, has been shown to have two 
regimes that can be approximated by two timescales, t\ 
and T2, which are about one order of magnitude apart, 
and are closely related to the r, t m , and the intercon- 
version rates: t\ = t(1 + 3/2cj), or t\ — > t m , and 
t-2 = t(1/2 + lo) + 1/ujq. From these relations the in- 
terconversion rates can be deduced when analyzing ex- 
perimental data. 
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APPENDIX A 

In this appendix we introduce the formal solution of 
Eq. 6 and define the symbols used in next derivations. 
In general, S(t) for a discrete system is given by summing 
the elements of the vector that solves the ME, 



S(t) = U 2n ~E,e ut ~E- 1 P (0|2n) 



(Al) 



Here U 2n is the summation row vector of 2n dimensions, 
P(0|2n) is the initial condition column vector, 



[P(0\2n)]j = {P A fiS x>j + PB,oS x+ n,j), 



(A2) 



where x is the initial state of the translocation process. 
The definite negative real part eigenvalues matrix, D, 
is obtained through the similarity transformation: D = 
E HE, where H is the matrix given on the right hand 
side of Eq. 6, and E and E _1 are the eigenvectors matrix, 
and its inverse, of H. 

APPENDIX B 

Here we calculate formally the MFPT < t >. The m 
moment of F(t) is given by: < r m > = J °° t m F(t)dt = 

m!f7 2 n(-H)- ro P(0|27i). To calculate the inverse of 
the propagation matrix H, which is given on the right 
hand side of Eq. 6 , we use the p r ojecti on operator 
i|Klafter and Silbev L 1198(1: IZwanzigt l200lft : QHQ = 

Hqq = A — U>A, Hq Z = Wb, H Z q = Wa, H Z z = 

B— ub, and the identity, I = HM, and obtain M blocks: 

Mqq = [A QQ - Aqz(Azz) _1 A Z q] _1 

= A" 1 C(B - w B ) 
M QZ = [A ZQ - A Z z(Aqz) _1 Aqq]- 1 

= -A^Cwb, (Bl) 

where Mzq and Mzz are obtained in a similar way. 
Now, we can write the m moment vector of F(t) as 



< r m > =m!(-M) m P(0|2n) 
where M is given by 



M 



A" 1 C(B-o;b) -A^Clob 
-A^Cua A- 1 C(A-w A ) 



(B2) 



(B3) 



and C=(B — wb — Aoja) For m 
obtain the MFPT vector: 



1 in Eq. B2 we 



< t >- 



( -A- x C{P A>0 ^-^)P{Q\n) 
\ -A- l C(P B , o A-cj A )P(0\n) 



(B4) 



where [P (0\n)]j=8 x j- Summing < r > elements by 
using the summation row vector of n dimensions U n , 
results in: 

< t >= -U n CP(Q\n)(XP A , + Pb,o) + 

+ U n A- 1 CP(0\n)(uj A +uj B ). (B5) 

Note that the MFPT of the single A conformation, r, is: 
t=— U nA^ 1 P (0|n), which has a similar form as of the 
first term in Eq. B5 when choosing C _1 as the propaga- 
tion matrix. 

It is easy to verify that for A=l, < r > reduces to the 
MFPT of the single conformation case, f . Rewriting Eq. 
B5 as 

< r >= -C^A^CIPa.oB + P B , A - 

-oj A - w B ] P(0|n), (B6) 

and substituting A=l we find that 

< t >= - UnA- 1 P(0\n) = t. (B7) 

APPENDIX C 

In order to obtain an explicit expression for the MFPT 
of the translocation in a weak field limit, we first rewrite 
Eq. B5 as 

< t >= ?(AP A ,o + Pb,o) + ^ 2 (^a + Wfl), (CI) 

where f = -U n CP(0\n), and ^ = J7 n A -1 CP(0|n). 
We can further rewrite r as r = t s , x ■ t s , x = 

— (C) SyX defines the mean residence time (MRT) spent 
in state s before exiting the channel, given that the pro - 
cess started at st ate x (|B_arzHMm_and Klafter L IT? 
and has the form llHuang and McColl 111997ft : 



'(C), 



A(/i s )A(/i 



— x 



) r x 



A(/i)A(^ 



s<x, (C2) 



where (C) S)K for s > x, is obtained when exchanging x 
for s and r for / in Eq. C2. Here h±=[\ ± (1 -4W) 1/2 ]/2, 
r=ap + , l—ap_, a~Xk/k, and k=Xk + + Aw^. Thus, 
we find that t s x is a function of the parameter a = [1 + 
{u) A + (jJs/A)/fc] _1 , which obeys < a < 1, and is a 
measure of the difference between r and r. Using a; w 1 
and u> A /k w 10^ 3 leads to a ~ 1 given V/Vc > 1, and 
accordingly to 

(C)., a = (A- 1 ),,^/^ - (A- 1 ). )(e /A. (C3) 
Eq. C3 implies 



r/A. 



(C4) 



To obtain an explicit expression for r, it is conve- 
nient to use the independence approximation and re- 
place Pj.j-i and kj by state independent terms: p+ = 
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[1 + e(~ v / Vc+1 ^]~ 1 and k. This approximation, which is 
valid for large polymers, N > d, and becomes more accu- 
rate as N increases, leads to a+ = p+fc, a _ = (1 — p+)k, 
so that ijFlomenbom and Klafter l l2003a|) : 

A{p n+1 - x )p x + x - A{p*)p n _ +1 - X (n +l-x) 
T ~ /cApA(p«+ 1 ) ' ( ' 

where A(p s ) = p s + — p s _. In the limit of a not too large 
field, V/Vc > 1, Eq. C5 reduces to 



2x£ p b 2 d> x 1 
z\q\{l + \/d)V -V c 



(C6) 



To compute a 2 we rewrite a 2 as c 2 = £" =1 t s t SiX , where 
t s is given by Eq. C5 for x = s, and r SjX is given by 

Eq. C2. For a w 1 we have cr 2 = r 2 /2A, where t 2 is 
the second moment of F(t) for the single ^4 conformation 
case. The calculation of r 2 /2 = J2™=i ~t s t s ,x yields in 
the weak field limit V/V c > 1: 



1 x(x-l) xy x (l - y) - y(l - y x ) 



2 (kAp) 2 1 2 
where y = p_/p + and 



-in, 

(C7) 



7J = 



?T X - 1 

y-n _ ! 



- 1 



n+1 — x 



1-^ 



+ n+ 1 
n + .t 



(C8) 



Noticing that 77 represents the non-translocation events 
and vanishes for V/Vc ^ 1 as y n ~ x , we rewrite Eq. C7 
up to a leading term in x as: 



— 1) 
2(kAp) 2 ' 



(C9) 



Using r ~ (x/kAp) valid for V/Vc > 1 
l|Flomenbom and Klafterl l2003a(l - Eq. C9 yields for a 
leading order in x 



—2 



/2A. 



(CIO) 



Substituting Eq. C4 and Eq. CIO into Eq. CI, Eq. 8 is 
obtained. 



which can be written as 

< t >= 7(1 + uj) + P b . /ujb. 



(D2) 



Note that experiments suggest that Ao <C 1, which leads 
to Aofc < lub, whereas Xok <C is used for simplifica- 
tion, and enables the substitution of A=0 in Eq. B5. 

To compute < t 2 >, we have to calculate the blocks of 
M 2 : 



M. 



QQ - = A 2 (1+uj); 
M QZ 2 = A~ 2 (l + u) - AT l /u B ; 

2 _ A -2, 



M zz 2 = A- 2 u(l + w) - A- 1 2uj/ujb + (D3) 

Substituting Eq. D3 into Eq. B2 and summing the vector 
elements, we obtain 



< T Z > 



— {1+uj) 2 + t(1 + 2lo) 
2 



^b.o , Pb.o 



(D4) 



To get the relaxation timescales of G ap (t), t\ and ra, we 
match < r > and < t 2 > obtained from 



F ap (t) 



Pa,o 
n 



t/Tl 



T 2 



t/r 2 



(D5) 



to the corresponding moments obtained from Eq. D2 and 
Eq. D4. This procedure yields: 



n 



<t > -Pb,qt 2 
Pa,o 



and 



t 2 =< r > 



Pa,o 



— <T> 



1/2 



(D6) 



(D7) 



Pb.o V 2 

Substituting Eq. D2 and Eq. D4 into Eq. D7 results in 
P Bfi , iVo 



= (l+u;) + 



1 — T- 



,2 



B 



Pa,o Pa,oPb,o 



■T 2 



1/2 



(D8) 



Expanding the square root in Eq. D8 to leading order 
and using Eq. D6, Eq. 11 is obtained. 
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